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Abstract 



We propose a new explanation for the longstanding puzzle of the tiny branch- 
ing fraction of ip' — > pir relative to that for J/tj} — * pn. In the case of J/ip, we 
argue that this decay is dominated by a higher Fock state in which the cc pair 
is in a color-octet 3 Si state and via annihilate process cc — > qq. In the case of 
the ip 1 ', we argue that the probability for the cc pair in this higher Fock state 
to be close enough to annihilate is suppressed by a dynamical effect related 
to the small energy gap between the mass of the ip' and the DD threshold. 
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Typeset using REVTgX 



A long-standing mystery of charmonium physics is the "p-7r puzzle" of J /if) and if)' 
decays. These particles are nonrelativistic bound states of a charm quark and its antiquark. 
Their decays into light hadrons are believed to be dominated by the annihilation of the cc 
pair into three gluons. In order to annihilate, the c and c must have a separation of order 
l/m c , which is much smaller than the size of the charmonium state. Thus the annihilation 
amplitude for an S-wave state like J /if) or if)' must be proportional to the wavef unction 
at the origin, if)(r = 0). The width for decay into any specific final state h consisting of 
light hadrons is therefore proportional to |?/>(0)| 2 . The width for decay into e + e~ is also 
proportional to |^(0)| 2 . This leads to the simple prediction that the ratio of the branching 
fractions for if)' and J /if) is given by the "15% rule": 

* s WJ^h = Q " = < 14 - 7±2 ' 3 > % - « 

The p-7r puzzle is that the prediction (1) is severely violated in the pre and several other decay 
channels. The first evidence for this effect was presented by the Mark II collaboration in 1983 
IJ. They found that Q pn < 0.6% and Qk*k < 2%. Recent data from the BES collaboration 
has made the puzzle even sharper. They obtained Q pn < 0.23%, Qk*+k- < 0.64%, and 
Qk*°k° — (1-7 ± 0.6)% 0. Thus the suppression of Q relative to the 15% rule is about 10 
for K*°K° and greater than 65 for pn. 

A summary of the proposed explanations of the p-7r puzzle has recently been given by 
Chao P . Hou and Soni ||] suggested that J /if) — > pn is enhanced by a mixing of the J / if) with 
a glueball O that decays to pn. Brodsky, Lepage, and Tuan || emphasized that J /if) — > pn 
violates the helicity selection rule of perturbative QCD, and argued that the data requires 
O to be narrow and nearly degenerate with the J/if>. Present data from BES constrains the 
mass and width of the glueball to the ranges \mo — mj/^l < 80 MeV and 4 MeV < To < 50 
MeV @. This mass is about 700 MeV lower than the lightest J PC = 1 glueball observed in 
lattice simulations of QCD without dynamical quarks J7|. There are explanations of the p-n 
puzzle that involve the dependence of the decay amplitude on the energy of the charmonium 
state. Karl and Roberts || suggested that the decay proceeds through cc — > qq followed 
by the fragmentation of the qq into pn. They argued that the fragmentation probability 
is an oscillatory function of the energy which could have a minimum near the mass of the 
if)' . Chaichian and Tornqvist pointed out that the suppression of ip' decays could be 
explained if the form factors for two-body decays fall exponentially with the energy as in 
the nonrelativistic quark model. There are other explanations of the p-7r puzzle that rely 
on the fact that there is a node in the radial wavefunction for if)', but not for J /if). Pinsky 
|I(J suggested that this node makes if)' — > pn a "hindered Ml transition" like J /if) — » ^7. 



Brodsky and Karliner |LT] suggested that the decay into p7i proceeds through intrinsic charm 
components of the p and it wavef unctions. They argued that the cc pair in the \udcc) Fock 
state of the p + or tt + has a nodeless radial wavefunction which gives it a larger overlap with 
J /if) than if)'. Finally, Li, Bugg, and Zou [0 have suggested that final-state interactions 
involving the rescattering of a\p and a,2p into pre could be important and might interfere 
destructively in the case of the if)'. 

In this Letter, we present a new explanation of the p-7r puzzle. We argue that the decay 
J /if> —>■ pir is dominated by a higher Fock state of the J / if) in which the cc is in a color-octet 
3 S\ state. The cc pair in this Fock state can annihilate via cc —>■ g* —>■ qq. The amplitude 
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for forming pir is dominated by the endpoint of the meson wavefunctions, where a single q 
or q carries most of the momentum of the meson. The suppression of if)' — > pre is attributed 
to a suppression of the cc wavefunction at the origin for the higher Fock state of the if)' . 
Such a suppression can arise from a dynamical effect associated with the small energy gap 
between the mass of the if)' and the DD threshold. 

It is convenient to analyze charmonium decays using nonrelativistic QCD (NRQCD) ||14|| , 
the effective field theory obtained by integrating out the energy scale of the charm quark 
mass m c . In Coulomb gauge, a charmonium state has a Fock state decomposition in which 
the probability of each Fock state scales as a definite power of v, the typical relative velocity 
of the charm quark. The dominant Fock state of the J /if) and the if)' is \cci( 3 Si)), whose 
probability V is of order 1. We denote the color state of the cc pair by a subscript (1 for 
color-singlet, 8 for color-octet) and we put the angular momentum quantum numbers in 
parentheses. The Fock states \cc 8 ( 3 Pj) + S), where S represents dynamical gluons or light 
quark- ant iquark pairs with energies of order m c v or less, have probability V ~ v 2 . The next 
most important Fock states include |cc 8 ( 1 S , ) + S) with V ~ v 3 (or ? ~ v 4 if perturbation 
theory is sufficiently accurate at the scale m c v) and |ccs( 3 Si) + S) with V ~ v 4 . 

The decay of the J / if) into light hadrons proceeds via the annihilation of the c and c, which 
can occur through any of the Fock states. This is expressed in the NRQCD factorization 
formula for the inclusive decay rate which includes the terms 

+ l ^4(O s ( 3 P a ))jn, + ^(O 8 (»S0>w + • • • • (2) 

where a s = a s (m c ). The matrix elements are expectation values in the J /if} of local gauge- 
invariant NRQCD operators that measure the inclusive probability of finding a cc in the 
J/ifj at the same point and in the color and angular-momentum state specified. The matrix 
element of the cci( 3 Si) term in (fj) is proportional to the square of the wavefunction at the 
origin and scales as v 3 . Its coefficient includes a term of order a 3 from cc — > ggg and a term 
of order a 2 from the electromagnetic annihilation process cc — > 7* — * qq. The color-octet 
terms in (§) represent contributions from higher Fock states. Their matrix elements scale 
like f 6 , v 7 and v 7 , respectively. Their coefficients are all of order a 2 and come from cc — > gg 
for the cc 8 ( 1 S'o) and cc 8 ( 3 -Po) terms and from cc — ^ g* — > gg for the cc 8 ( 3 S'i) term. Note that 
the coefficients of the color-octet matrix elements are two orders of magnitude larger than 
that of (Oi( 3 Si)) ju, which suggests that the higher Fock states may play a more important 
role in annihilation decays than is commonly believed. 

Bolz, Kroll, and Schuler have emphasized that contributions from higher Fock states 



should also be important in exclusive decays of charmonium into light hadrons |R|. They 
argued that the amplitude for a two-body annihilation decay such as J /if) — > pn satisfies a 
factorization formula which separates the scale m c from the lower momentum scales. The 
decay amplitude is expressed in terms of hard-scattering factors T that involve only the 
scale m c , initial-state factors I that involve scales of order m c v and lower, and final-state 
factors T that involve only the scale Aqcd- If m c was asymptotically large, the dominant 
terms in the factorization formula would have the minimal number of partons involved in the 
hard scattering. Terms involving additional soft partons in the initial state are suppressed 



3 



by powers of v. Terms involving additional hard partons in the final state are suppressed 
by powers of Aqcd/^c- Thus the asymptotic factorization formula has the schematic form 

T(J/ip -> P+7T-) = £ [f (cc^Sx) -> ud + du)F p+ (ud) ^-(du)] Tj / ^{cc l { z S l )). (3) 

ud,du 

The initial-state factor is proportional to the wavefunction at the origin for the \cci( 3 Si)) 
state. The lowest order contributions to the hard-scattering factor are of order a 3 from 
the annihilation of cc into three virtual gluons and aa s from the annihilation into a virtual 



photon. As pointed out by Brodsky and Lepage ||15|| , the leading term in the asymptotic fac- 
torization formula (0) is strongly constrained by the vector character of the QCD interaction 
between quarks and gluons. It vanishes unless the sum of the helicities of the mesons is zero. 
Rotational symmetry then requires the angular distribution to be 1 — cos 2 9. This helicity 
selection rule is violated by the decay into pit. Parity and rotational symmetry require the 
helicity of the p to be ±1 and the angular distribution to be 1 + cos 2 9. Since the helicity 
of the pion is 0, the helicity selection rule is violated. Thus the amplitude for J/ip —>■ pir is 
suppressed by a factor of Aqcd/^c relative to that for some other mesons. 

Since the charm quark mass m c is less than an order of magnitude larger than Aqcd, 
there can be large corrections to the asymptotic decay amplitude (|3]). In particular, there 
can be significant regions of phase space in which some of the gluons involved in the hard 
scattering are relatively soft. It might therefore be more appropriate to absorb them into 
the initial-state or final-state factors. In this case, not all of the soft partons in S need be 
involved in the hard scattering, and not all the partons that form the p + and ir~ need be 
produced by the hard scattering. For example, there can be a contribution from the Fock 
state |cc 8 ( 3 Si) + S) that involves the hard-scattering process cc — > uu. This produces a state 
\uu + S) that consists of the soft partons S together with a u and a u that are back-to-back 
and whose momenta are approximately m c . Such a state has a nonzero overlap with the final 
state |p + 7r~). The overlap comes from the endpoints of the meson wavef unctions, in which 
most of the momentum of the p + is carried by the u and most of the momentum of the 7r~ is 
carried by the u. There is also an endpoint contribution involving the dominant \cci( 3 Si)) 
Fock state, which can annihilate into qq through a virtual photon. These contributions to 
the T-matrix element can be written schematically in the form 

r(J/V> - P+7T-) = £ f(cc 8 ( 3 5 1 ) - qq) £ [F^-(qq + S) J J/ ^(cc 8 ( 3 S 1 ) + S) 
qq s 

+ £ ?(cc x <?S x ) -> qq) F p+ „-{qq) l J/i) {cc l ^S l )) , (4) 
qq 

where the sum over qq includes uu and dd. In the factorization formula (0) for inclusive 
decays, the two contributions on the right side of (|j) contribute to the cc 8 ( 3 Si) and cci( 3 Si) 
terms, respectively. 

The endpoint contribution in (^) leads to a definite angular distribution. Since the q 
and q carry most of the momenta of the mesons, the angular distribution of the mesons will 
follow that of the q and q, which is 1 + cos 2 9. Thus @ will contribute most strongly to form 
factors which allow the angular distribution 1 + cos 2 9. It will also contribute most strongly 
to decays into mesons like p and ir for which most of the momentum can be carried by a 
single q or q. There are also endpoint contributions involving the hard-scattering process 
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ccs( 3 Pj) — > gg, which produces a pair of hard gluons with the angular distribution 2 + cos 2 9, 
and ccg( 1 5'o) — > gg, for which the angular distribution is isotropic. Their contributions to 
J/i/j — > pix are suppressed by the small probabilities for most of the momentum of the p or 
7i to be carried by a single gluon and by the mismatch between the angular distribution of 
the gluons and the 1 + cos 2 9 distribution of the pn. 

We argue that the color-octet term in (f|) may actually dominate the decay rate for 
J/ip — ► p7r. We compare the various factors in that term with those in the asymptotic 
expression (|3|). The hard-scattering factor T in the color-octet term in (f|) is only of order 
a s , compared to a 3 in (|3|). The suppression from the initial-state factor in including 
the sum over S, might be as little as a factor of v 2 relative to that in (0) if the additional 
soft partons S consist of a single gluon with momentum of order m c v 2 . As for the final- 
state factors, (Q) is suppressed by the endpoints of the meson wavef unctions, while (|3|) is 
suppressed by Aq CD /m c from the violation of the helicity selection rule. Considering the 
various suppression factors, it is certainly plausible that the cc 8 ( 3 Si) term in ([|) could 
dominate over the leading contribution from the cci( 3 Si) Fock state. 

We have argued that the decay J/ip —>■ pn may be dominated by the annihilation of the 
cc pair in the \cc$( 3 Si) + S) Fock state via cc — > qq. If this is true, then the p-n puzzle 
can be explained by a suppression of this decay mechanism in the case of the ip' . This 
suppression can arise from the initial-state factor Ij/^p(cc s ( 3 Si) + S) if the wavefunction for 
the \cc$( 3 Si) + S) Fock state is suppressed in the region in which the separation of the cc 
is less than or of order l/m c . Note that it does not require a suppression of the probability 
for the higher Fock state, but just a shift in the probability away from the region in which 
the c and c are close enough to annihilate. A possible mechanism for this suppression is a 
dynamical effect related to the small energy gap between the mass of the ip' and the DD 
threshold. In the Born-Oppenheimer approximation |16|| , the cc pair in the dominant \cc\) 
Fock state moves adiabatically in response to a potential V\(R) given by the minimal energy 
of QCD in the presence of a color-singlet cc pair with fixed separation R. Similarly, the cc 
pair in the \cc$( 3 Si) + S) Fock state moves adiabatically in response to a potential Vs(R) 
given by the minimal energy of the soft modes S in the presence of a color-octet cc pair 
with fixed separation R. At short distances, this potential approaches a repulsive Coulomb 
potential a s /(QR). At long distances, the minimal energy state consists of D and D mesons 
separated by a distance R, and V$(R) therefore approches a constant 2(Md — m c ) equal to 
the energy of the DD threshold. A charmonium state spends most of the time on the color- 
singlet adiabatic surface, but it occasionally makes a transition to the color-octet adiabatic 
surface. Since the J/ip is 640 MeV below the DD threshold, a cc pair on the color-octet 
adiabatic surface is far off the energy shell. The time spent by the J/if> on this surface is 
therefore too short for the cc pair to respond to the repulsive short- distance potential. The 
shape of the cc wavefunction for the |ccs + S) Fock state should therefore be similar to that 
of the |cc"x) Fock state, which peaks at the origin. However the mass of the ip' is only 43 
MeV below D + D~ threshold and 53 MeV below D°D° threshold. A cc pair on the color- 
octet adiabatic surface can be very close to the energy shell. The ip' can therefore spend a 
sufficiently long time on this surface for the cc pair to respond to the repulsive short- distance 
potential. This response can lead to a significant suppression of the cc wavefunction at the 
origin for the |ccg + S). 

If the initial-state factor Xj/^,(cc 8 ( 3 S'i) + S) in (|J) is suppressed for all soft partons S, the 
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suppression can be expressed in the form of a relation between the NRQCD matrix elements 



in 



This inequality can be tested by calculating the matrix elements using Monte Carlo simula- 
tions of lattice NRQCD. Since the if)' is so close to the DD threshold, it would be essential 
to include dynamical light quarks in the simulations. 

Our proposal leads to a prediction for the flavor- dependence of the suppression of the 



decays of if)' into vector /pseudoscalar final states. Bramon, Escribano, and Scadron [17 
have analyzed the decays J /if) — > VP assuming that the decay amplitude is the sum of a 
flavor-connected amplitude g, a flavor-disconnected amplitude rg, and an isospin-violating 
amplitude e. They allowed for violations of SU(3) flavor symmetry through parameters s 
and x — 1. Two of these amplitudes are 

T(J/ip -> pn) = (g + e) e j/ip x e p ■ p w , (6) 
T(J/if) - K*°K°) = [(1 - s)g - (1 + x)e\ e w x e K * ■ p R . (7) 



The amplitudes for the other VP decay modes are given in Ref. [17]]. The authors give two 
sets of parameters that fit the existing data, one with x = 1 and the other with x = 0.64. 
Both sets have e comparable in magnitude to rg and an order of magnitude smaller than 
g. If the decay J /if) — > pi: is dominated by endpoint contributions, we can identify g and e 
with the two terms in (^j). While rg may also have endpoint contributions from cc — > gg, we 
assume for simplicity that it is dominated by subasymptotic contributions from the cci( 3 5i) 
Fock state. The amplitudes for the decays if/ — > VP can be expressed in a similar way in 
terms of amplitudes g', e', and (rg)'. Our explanation of the p-7r puzzle implies that \g'\ is 
much smaller than \g\, and that e' and (rg)' differ from e and rg only by the factor required 
by the 15% rule. The unknown amplitude g' is constrained by the BES data on if)' — ► pn 
and if)' — > K*°K°. The upper bound on B(if>' pix) gives an upper bound on \g' + e'| 2 , 
which implies that g' lies in a circle in the complex (/-plane. The BES measurement of 
B(ip' -> gives an allowed range for |(1 — s)g' — (1 + x)e'| 2 , which constrains to 

an annulus. The intersection of the interior of the circle with the annulus is the allowed 
region for g' . By varying g' over that region and taking into account the uncertainties in the 



parameters of Ref. ||17|| , we obtain the predictions for Qyp in Table 1. Measurements of the 
if)' branching fractions consistent with these predictions would imply that the suppression 
of the vector /pseudoscalar decays is due to the suppression of g' . This would lend support 
to our explanation of the p-n puzzle. 

Our proposal also has implications for the angular distributions of other two-body decay 
modes. In general, the angular distribution must have the form 1 + acos 2 #, with —1 < 
a < +1. Our solution to the p-ir puzzle is based on the suppression of a contribution to if)' 
decays that gives the angular distribution 1 + cos 2 9. Thus the parameter a for any two-body 
decay of the if)' should be less than or equal to a for the corresponding J /if) decay. 

A solution to the p-n puzzle should also be able to explain the pattern of suppression 
for various J PC states with the same flavors. A preliminary measurement of the axial- 
vector/pseudoscalar decay mode if)' — > biir by the BES collaboration gives Qbm — (24 ± 



6 



7)%, consistent with no suppression relative to the 15% rule. A preliminary measurement 



of the vector/tensor decay mode if)' — > pa 2 [18| gives Q pa2 = (2.9 ± 1.6)%, which, though 
suppressed relative to the 15% rule, is an order of magnitude larger than Q p7T . This pattern 
can be explained by also taking into account the orbital-angular-momentum selection rule 
for exclusive amplitudes in perturbative QCD [[0J. The decay modes biir and pa 2 both have 
form factors that are allowed by the helicity selection rule. They also both have form factors 
that violate the helicity selection rule, but are compatible with an endpoint contribution 
from cc — > qq. However, in the case of bin, the endpoint contribution is further suppressed 
by the violation of the orbital-angular-momentum selection rule. A thorough analysis of J/ip 
and ip' decays into axial- vector /pseudoscalar and vector/tensor final states will be presented 
elsewhere. 

In conclusion, we have proposed a new explanation of the p-n puzzle. We suggest that 
the decay J/ip — > pit is dominated by a Fock state in which the cc is in a color-octet 
3 S\ state which decays via cc — > qq. The suppression of this decay mode for the ip' is 
attributed to a dynamical effect that suppresses the cc wavefunction at the origin for Fock 
states that contain a color-octet cc pair. Our explanation for the p-7r puzzle can be tested 
by studying the flavor dependence of the two-body decay modes of the J/ip and ip', their 
angular distributions, and their dependence on the J PC quantum numbers of the final-state 
mesons. 
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TABLES 



T/n 

V r 


X 


= 1 


x - 


- 0.64 


pir 


- 


0.25 


- 


0.25 


K*°K° + c.c. 


1.2 


- 2.0 


1.2 


- 3.0 


K* + K~ + c.c. 


- 


0.36 


- 


0.52 


LOT] 


- 


- 1.6 


- 


- 1.6 


LOT]' 


10 


- 51 


12 


- 55 


4>T] 


0.8 


- 3.6 


0.4 


- 3.0 


4>T]' 


0.7 


- 2.5 


0.5 


- 2.2 


prj 


14 


- 22 


14 


- 22 


pv' 


12 


- 20 


13 


- 21 


LVTT 


11 


- 17 


11 


- 17 



Table 1. Predictions for Qyp in units of 1% for all the vector /pseudoscalar final states. The 
values for pir and K*°K° + c.c. were used as input. The columns labelled x = 1 and x = 0.64 
correspond to the two parameter sets of Ref. 
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